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Abstract 

□ 

In Aldous and Shields (1988), a model for a rooted, growing random binary tree was presented. 
For some c > 0, an external vertex splits at rate c~' (and becomes internal) if its distance 
from the root (depth) is i. For c > 1, we reanalyse the tree profile, i.e. the numbers of external 
vertices in depth i = 1,2, .... Our main result are concrete formulas for the expectation and 
covariance-structure of the profile. In addition, we present the application of the model to 
cellular ageing. Here, we assume that nodes in depth h + 1 are senescent, i.e. do not split. 
We obtain a limit result for the proportion of non-senescent vertices for large h. 

1 Introduction 

Trees arise in several applied sciences: In linguistics and biology, trees describe the relationship 
of items (languages, species) and in computer science, trees are used as data structures, e.g. 
for sorting. Randomizing the input leads to random trees, which are object of a large body of 
research. For applications in biology, see e.g. [5J[T3]. Here, important examples are trees arising 
from branching processes (e.g. Yule trees). In computer science, prominent examples are search 
trees; see e.g. [T5J[TU]. 

In this note, we are concerned with an application of random trees in cellular biology. In 
the 1960s it was known that eukaryotic cells have a limited replication capacity ([Hj). The 
number of generations until cells do not proliferate any more is today known as the Hayflick 
limit and the phenomenon that cells loose their ability to proliferate is called cellular senescence. 
The molecular basis for cellular senescence were uncovered starting in the 1970s. A theory was 
developed which argued that during each round of replication, the telomeres (which are the end 
part of each chromosome) are shortened due to physical constraints of the DNA copying mechanism 
(|20j). In humans, these telomeres are a multiple (i.e. more than lOOG-fold) repetition of the base 
pairs TTAGGG and up to 200 bases are lost in each replication round ([TTI). Most importantly, 
telomeres have a stabilizing effect on the DNA. The DNA repair mechanism of a cell must be able 
to distinguish between usual DNA breaks (which it is assumed to repair) and the telomeres (which 
it is assumed to ignore). Hence, when telomeres become shorter this stabilizing effect seizes and 
ageing occurs. It can be observed that telomeres shrink from 15 kilobases at birth to less than 5 
kilobases during a lifetime ( [24 ) . However, the enzyme telomerase is known to be able to decrease 
the loss of telomeres during replication. This enzyme has been found to be active in stem cells 
and cancer cells, which both are cell types with an (almost) unbounded replication potential. The 
deeper understanding of the role of telomeres and telomerase is an active field of research because 
of the medical implications for ageing and cancer. In particular, it was awarded the Nobel prize 
in medicine in 2009 (US]). 
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We study the model of random trees introduced in Aldous and Shields in jTj (hereafter referred 
to as [AS]) and extend it for an application to cellular ageing. Given some c > and a full binary 
tree T, the model introduced in [AS] describes the evolution of the vertices of the tree. Here, 
we distinguish internal, external and prospective vertices. At t = 0, the root is the only external 
vertex (and there are no internal vertices). An external vertex m G T in depth \u\ becomes internal 
at rate c"'"'. At the time it becomes internal, the two daughter vertices in depth |u| + 1 become 
external. We present our result on the profile of the Aldous-Shields model in Theorem [T] 

For our application to cellular senescence, we will analyze a relative of the Aldous-Shields 
model for c > 1. Here, a critical depth h is fixed, and only external vertices in depth at most 
h can become internal. External vertices in depth h + I never become internal. Here, external 
vertices can be thought of as cells. The depth of a vertex is the number of generations from the 
first cell. Vertices in depth at most h represent proliferating cells, because they are able to produce 
offspring (i.e. daughter cells). Vertices in depth h + 1 represent senescent cells. This model has 
two features, which appear to be realistic in cellular senescence. First, the rate of cell proliferation 
decreases with the generation of a cell, parameterized by c > 1. Second, cells which have already 
split too often loose their ability to proliferate at all. For this model, we obtain a limit result for 
the frequency of proliferating cells in Theorem [2j 

The paper is organized as follows: In Section [2j we state our results on the Aldous-Shields 
model. The application to cellular senescence is carried out in Section [3j where we also give an 
overview of other models for cellular senescence in the literature. Section |4] contains the proofs 
for our results on the Aldous-Shields model (Theorem [ij , and in Section [5j we give proofs for the 
results on the model of cellular ageing (Theorem [2]) . 

2 Model and results 

We start by introducing some notation. Let T be the complete binary tree, given through 

oo 
n=0 

and 

To = {0}, T„ = {0,1}" forn=l,2,... 

We refer to elements in T by vertices and identify u G T„ by a word of length n over the alphabet 
{0, 1}, whose ith letter is Ui, n > 1. The vertex is the root of the tree and vertex u €T has two 
daughter vertices, uO and ul. (We make the convention that 00 0,01 := 1.) For u E T we set 
|u| = n iff u G T„. 

We say that u is an ancestor of v if |m| < \v\ and there are ii, i|^|_|„| G {0,1} with v — 
uii ■ ■ -i^^^-^u]- The ancestor induces a transitive order relation in T, and we write m ^ u iff m is 
ancestor of v. 

Definition 2.1 (Aldous-Shields model). Fixe > 0. The (time- continuous) Aldous-Shields model 
with parameter c is a Markov jump process y = {Y{t))t>o, Y{t) ~ (i^ii(t))tteT with state space 
{0, 1}'''', starting in Y{0) — {lu=(ti)uei: ■ Given Y{t) = y E {0, 1}'"' and u E T with y„ = 1, it jumps 
to iyv)veT, given by 

{0, V = u, 
1, V = uO or V — ul, 
yv, else, 

at rate c~'"'. In this case, we say that vertex u splits. 
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Figure 1: Two realizations of the Aldous-Shields model at the time when 500 vertices are external 
for c = 1.05 (A) and c — 3 (B). Only depths 0, ...,10 are drawn and the external vertices are 
marked. 



Remark 2.2 (Internal and external vertices). Let y = {Y{t))t>o be the Aldous-Shields model 
and Y = Y{t) for some t > 0. It is important to note that the dynamics is such that any path 
0, ii, iii2, ■■■ G T with ii, Z2, ... £ {0, 1}, starting at the root, has exactly one element u with Y^ = 1. 
In particular, the sets 

{u : 3v : u ^ v,Yy — 1}, {u : Y^ = 1}, {u : 3v : v ^ u^Yy = 1} 

of internal, external and prospective vertices are disjunct. 

Definition 2.3 (Profile). Let y = {Y{t))t>o be the Aldous-Shields model and Y = Y{t) for some 
t>0. We define 

Xn := Yu, Xn :— 2 "A"„, (2-1) 
tieT„ 

the total number of external vertices and the relative proportion of external vertices in depth n, 
respectively. The vector {Xn)n=o,i,2,... is also called the profile and 

oo 

x = Y,x^ 

is the total number of external vertices. 

Remark 2.4 (Dependence on c). The behaviour of the Aldous-Shields model strongly depends on 
c. A larger c implies that the profile is more concentrated around certain depths. This is because 
a larger c means that external vertices in smaller depth have a higher chance to be the next to 
split. See Figure [T] for an illustration. 

Two values of c are of particular important in applications from computer science: for c = 1, 
and ii X = n, the set of external vertices is a binary search tree with n external vertices. For 
= 2 and X = n, the set of external vertices is a digital search tree with n external vertices; see 
e.g. [10]. 

Remark 2.5 (Relative frequencies). We observe that 



oo 
n=0 



(2.2) 
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for all t, almost surely. To see this, note that Xo{0) = 1 and X„(0) = for n > 0, i.e. (2.2) holds 
at i = 0. Additionally, assume that some u with \u\ — n splits at time t. Then, we have that 
X„(<) - X^{t-) = -2-" and Xn+i{t) - X„+i(t-) = 2 • 2-("+i) = 2"". In particular, every split 
leaves X^i^Lo^n unchanged which shows (2.2). 

Remark 2.6 (Notation). In our results, we will give asymptotics of moments of X„+j(ic") for 
large n. Generally, for two sequences {xn)n=i,2.... and {yn)n=i,2,..., which may depend on other 
parameters, we write 



n— f oo 



hm ^ = 1. 



Theorem 1 (Moments of the profile and their limits). Let c > 1. Define for k e 



au = {-ly 



(c-l)---(cfe-l)' 



■ ■ - c^ 



(c - 1) • • • (cfe - 1) 



(2.3) 



and ao = b{) = 1 (with the convention that an empty product equals 1). For negative k, we set 
ttk — bk — 0. Moreover, 



boo :=(n( 

1=1 

Then, for t > and i > —n, 

n 

E[X„+,(te")] -2"+^.^afe6„_fce 
Moreover, for i, i' G Z, i < i' , 



e (1;(X)). 



i + k. n— >-oo 



boo ■ 2"+' • J2 ^^e" 



(2.4) 



k=0 



2 \ 2n+i+i' 



2 - c2 V c 



[Xn+^{tc"),Xn+Atc")] « a^^,'•< 



2"r^^/2' 



c<V2 



(2.5) 



l2(c2-l)(c2-2) 



where 



b"^ ^ OkOk'e' 

k,k'=a 



tk+k' 



(2.6) 



Remark 2.7 (Convergence and covariance). 1. It is immediate from the Theorem that 

oo 

lim X„+,(te") = 600 • y flfee-^""^'* (2.7) 



fc=0 



in probability, for alH > and all c > 1. 

2. The covariances given in (2.5) show a phase transition at c = ^/2. Such a phase transition 
is already known from results by [AS] and [9]. However, these papers do not give explicit 
formulas for the covariance structure. 



3. Using (4.4), (4.6) and (4.10) it is also possible to obtain exact results for the covariance on 
the left hand side of (2.5). 
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Remark 2.8 (Connection to results from Aldous and Shields (1988)). In [AS], the evolution of 
the vector {X„{t))n=o,i,2,... is studied. In their Theorems, a law of large numbers if Xn{t) to some 
deterministic limit x„(t) is stated and proved using martingale methods. Their result implies that 



(2.7) holds even almost surely on compact time intervals. In particular, they claim that the limit 



Xi{t) of Xn+i{td^) must satisfy Xi{t) = Xi+i{ct), which clearly holds for the right hand side of 



(2.7). In addition, they show that a suitably rescaled process, 2"/^(X„+i(t) — Xi{t))t>o, converges 



as n — >■ oo weakly to a diffusion for c > \/2. The rescaling factor 2"/^ can also be seen from the 



above Theorem. Moreover, (2.5 1 shows that a convergence of c^{Xn+i{t) — Xi{t))t>o to a diffusion 



can be conjectured for 1 < c < v^. 

Remark 2.9 (Connection to work of Dean and Majumdar (2006)). In [9j, the total number of 
external vertices, X, was studied in the context of the Aldous-Shields model on an m-ary tree. In 
the binary case, a functional equation (their equation (2)) for the Laplace transform of X{t) was 
shown to hold true. This equation uses the following fact: Given that T is the random time of 
the first split in the model, it is clear that T is mean one exponential and, in addition, 

X{t) ^ lT>t + 1t<.(X'(*^) + ^"(^)) 

where X' and X" are independent of T and of each other and distributed like X. From their 
identity on Laplace transforms, [9j show the phase transition for the variance of the number of 
occupied vertices at c = v^, which is also seen from Theorem [I] 



3 Application: cellular ageing 

The first mathematical model for cellular senescence was given in [17 . It takes several biological 
facts into account. When DNA is copied, the double helix is unfolded and both strands of DNA 
are copied. Only in one of the two strands there are physical constraints by which the end of 
a chromosome cannot be perfectly copied. This shortening of telomeres is independent for all 
chromosomes. In [T7], a fixed length for telomeres which decreases by a fixed amount at each 
proliferation event for one of the daughter cells and proliferation occurs along a full binary tree 
is assumed. If the length of a telomere of one chromosome falls below a threshold, a cell cannot 
replicate any more and becomes senescent. This threshold takes the Hayflick limit into account, 
which states that a cell line can only life for a limited number of generations before it becomes 
senescent. 

The model by [T7] was extended in several directions. A stochastic amount of loss of telomeres 
was studied in [2]. In [3] and [19], the binary tree of proliferating cells from [17] was replaced 
by a branching model. In particular, [TH] took cell death into account, with different death rates 
above and below a critical threshold of telomere length. Age structure of cells (i.e. structure which 
phase of the cell cycle) is taken into account by [TTJ[T2]. Moreover, [J] extend the model of [T7] by 
explicitly taking telomerase activity (which is present in stem cells and cancer cells) into account. 

The idea to use the Aldous-Shields model for cellular ageing was influenced by the following 
recent results: 

1. In [7], a model is proposed which distinguishes two states of telomeres: capped and uncapped. 
Only in the capped state, proliferation of the cell is still possible. In somatic cells, an 
uncapped telomere cannot be transformed to the capped state any more leading to senescent 
cells; see the model of [22]. In stem and tumor cells, telomerase is (among other things) 
responsible for transitions from uncapped back to capped telomeres. Following [23], the 
transition rate of the uncapped to the capped state in stem cell decreases with shorter 
telomeres. 

2. In data, it has been observed that proliferating cells can behave differently. Motivated by 
data from [SlIHUH], it is argued in [21_ that the rate of proliferation decreases for shorter 
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telomeres. Their model produces a Gompertzian growth model which is known to fit to 
empirical data for somatic and tumor cells. 

In stem cells, the decreasing rate for an uncapped telomere to reeneter the capped state for shorter 
telomeres from [53] shows exactly the behaviour of the Aldous-Shields model: cells with a long 
replicative history proliferate slower. While [5T] use a linear decrease in replication rate, depending 
on telomere length, the Aldous-Shields model uses a geometric decay of the proliferation rate. 

Note that short telomeres can be seen as a form of damage. In [13] , models for cellular damage 
were introduced. In their model, cells inherit damage to the daughter cells. This model, as well 
as the Aldous-Shields model are among the analytically tractable ones. 

We state our model of cellular ageing: 

Definition 3.1. Fix h e N,r > and c> I and let T'' := Uj^toT,,. The process Z = {Z{t))t>o, 
where Z{t) = {Zu{t))ueTh is a Markov jump process with state space {0, iy\ starting in Z{Q) = 
(lu=0)tieT'' • Given Z(t) = z G {0, 1}''' and u G T'' \Tft^i with z„ = 1, it jumps to {zv)v(zj, given 
by 

0, V ^ u, 

1, V ^ uO or V ^ ul, 

at rate re"'"'. (Note that vertices u G Th+i do not split.) 

Informally, every external vertex u in this process represents a cell. If |m| — n, we say that 
the cell is in generation n. The process starts with a single mother cell. It proliferates at rate r. 
All cells up to generation h from the mother cell follow the usual dynamics of the Aldous-Shields 
model (with time rescaled by a factor of r), such that cells in generation n proliferate at rate 
r ■ c~". If a cell is in generation h + 1 from the mother cell, its telomeres have reached the Hayflick 
limit and the cell is not able to proliferate any more. 

Definition 3.2 (Relative frequency of proliferating cells). In applications, the relative frequency 
of proliferating cells, 

Lit) ■= , , (3.1) 

^' zp{t) + zHt) 



with 



+1 



is of particular importance. Here, Z^it) and Z^{t) is the number of proliferating and senescent 
cells at time t, respectively. 



Theorem 2 (Frequency of replicating cells). For Z and L as in Definition 3.1 and 3.2 



lim Lite'') = ^oo ^oT'"" . , , (3.2) 

in probability, for all t > Q. 

Remark 3.3 (Simulations). In our model for cellular senescence, Theorem[2]describes the decrease 
in the frequency of proliferating cells. This frequency has been measured empirically; see e.g. [H 
Figure 5] and [3| Figure 2]. As can be seen from the Theorem, every c gives a specific curve of 
decrease; see also Figure [2] These curves can be fit to data in order to estimate c. As the figure 
shows, the limiting result of Theorem [2] already gives a good fit for simulations which use h — 20. 
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o 



o 



o 
o 




Figure 2: The decrease of the frequency of prohferating cells strongly depends on c. The figure 



shows simulations for c = 1.2, 1.3, 1.5 and h = 20. The grey lines show the limit result (3.2) from 
Theorem 2 with r = 1. 



4 Proof of Theorem [T] 
4.1 Preliminaries 

The key ingredient in the proof is the quantity 



yn{t) := E[ro„(i)] with 0„ 0_^. 

n times 



(4.1) 



Note that the dynamics of (i1]„)n=o,i,2,... (taking values in {0, ijlO'i^^,...} ^j^j^ Qjjjy Qj^g entry being 
1) is autonomous. It is given by the following rule: if Yo^{t) = 1, then, at rate c~" a transition 
occurs to the configuration Yo^{t) = 0, yo„^i(i) = 1. From the dynamics of the Aldous-Shields 
model, it is clear that y(t) = {yn{t))n=o,i,2.... follows the differential equations 



with y_i — 0; compare (2.1) in [AS]. We rewrite the equation to obtain 

y ^ 



(4.2) 



with 



/-I 
1 



A = 



— c 
^-1 



— c 

^-4 



V 



■7 



Our first Lemma provides essential facts about the matrix A. Recall Ofc and bk from (2.3) 
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Lemma 4.1. Let E — (eij)ij=o,i,2,... o-nd F — (/ij)ij=o,i,2,... be given by 

^ij — ^i—j 1 fij — bi—j ■ 

The matrix E contains the eigenvectors of A to the eigenvalues Xi — — c^', i — 0,1, 2, ... and E 
and F are inverse to each other. 

Proof. To see that E contains the eigenvectors of A, note that for i > j 

{AE)ij = c"'+^a,_i_j - c^''a^-J = (c"'+^^^ — ^ c"') a,^j = -c^^Oi^j. 

To see that E and F are inverse to each other, it follows from the definition of at and bk that 

{FE)u = ao ■ bo = 1 and {FE)ij = for i < j. 
For i > J, we set n := i — j > and obtain 

i n+j n 

{FE)i,j — ^ bi^k^k-j — ^ bj^n^kak-j — ^ bn-kO-k 

k—j k—j k—0 

^ ^ . . . k ^k 

^ (c-l)---(c«-'=-l)-(c-l)---(c'=-l) 

" ^l + --- + (ji-/£-l) 

= (c - 1) . . . (c« - 1) (c - 1) • • • (c"-''^ - 1) ^"'^^ - 1) ■ ■ ■ (c" - 1) 

= ( u 1, E(-l)'^^f^7^^(^""'^' - 1) • • • (C" - 1), 

(c- 1)- •• (c« - 1) ^ (c-l)---(c*^-l) 

where we have reversed the order of the summands in the last equality. We rewrite (c' — 1) • • • (c' 
I' 

I) = Yl (tf - 1) = 1 for ; > /' and claim that for aU n > 

j=i 

E(-l)'7 .. (^"'"^^ - 1) • • • (c" - 1) = (4.3) 

^ (c-l)---(c'^-l) 

which implies that F and i? are inverse to each other. We use induction and note that the assertion 
is clear for n = 1. Given it is true for n, we have 

n+l (J) 

Yi-^f-, ^-^i. r(c""''^+' - 1) • • • (c"+i - 1) 



fc=i -LJ-'-lc ij 

n+l 

= E(-l)'7 u\k l)(c"-'=+2 _ 1) ... (c" - 1) 

^ (c-l)---(c'^-l) 

n+l C""^) 

= E(-l)(-l)'-' (e _ 1) . . . (c^-i 1 i)(e^ „ 1) i'^' - - 1) ... (c" - 1) 

n+l ^=-1) 

= E(-l)'' \,_i).'..(,.-i_i) (^"'''"^^' - 1) ... (c- - 1) = 

where we have used the induction hypothesis in the second and in the last equality. Hence, we 
have shown (4.3) and the proof is complete. □ 
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4.2 First order structure; proof of (2.4) 



By linearity, we can now explicitly solve (14. 21) using Lemma 4.1 Since E contains the eigenvalues 



of A and F is inverse to E, we immediately write, using y{t) :~ (?/o(^), yi(^), J/2(i), •■•) and D := 
diag(-l,-c-i,-c-2,...) 

n 



k=0 



Ti=0,l,2, 



because z(0) = (1, 0, 0, ...) as the process starts with {Y{0)) = (l„=0)„gT- We have shown the first 



part of (2.4 1 and in order to prove the second part, fix i, k, t and note that n i— )■ \ak\hn+i-kG 



is increasing with a summable limit. Hence, by dominated convergence. 



y„+i(te") = akbn+i-ke 



k=0 



k=0 



4.3 Second order structure; proof of (2.5) 



Now we come to the second order structure. Similar to the definition of y„ in (4.f ), we set for 
n < n' 

Vn.n'.mit) COW[Yo„it),Yo^i^,_Jt)] with 0„1„'_„ 0^ . 

m times n'—m times 

In order to see the connection of yn,n',m{t) and COV[X„(t),X„/(t)], we define the depth of the 
most recent common ancestor of u, w' as 

Muy := sup{|'y| : w ^ u, w ^ w'}, u, v! G T, 

where v ^ u if v ^ u or v = u. Let U, U' be two random variables, where U is uniformly distributed 
on Tn and U' uniformly distributed on T„/, independent of all the rest. The distribution of Mu^jjt 
is given by (recall n < n') 

¥[Mu,u' =m\= 2-«"'+i)^"), m = 0,...,n. 

{t),x,At)]^ E inyu{t)Yu'{t)]-E[Yu{t)]-E[Yu>{t)]) 

= 2''+'''COY[Yu{t),Yu'{t)] 

= 2"+"'(E[COV[yt/(t),n/'(i)|Ma,a']] + COW[E[Yu\Mu,u'],nYu'\Mu,u']\) (4.4) 
Yuit),Yu'it)\Mu,u']] 



We write 



= 2"+"'E 



= 2"+"' 2-«'"+i)^")y„,„,,„(t). 



7n— 

The second to last equahty holds as E[y[/|A%^t/'] = E[Yu] and E[Y;7/ |Af[/^[//] = E[Yt//]. In order to 
use the last expression, note that for yn,n',m the vertices 0„ and Om^n'-m have the vertex 0^ as 
their most recent common ancestor and so Afo„,Omi„/_„ = iti- We let be the last time Yq^ is 
external, respectively the time when Yq^^-^ becomes external, i.e. Tm is the sum of exponentials 
with parameters 1, c~^, c~^, c~™. Hence, 



El 



n 



+A 

1=1 1=1 



Ac'' 



(4.5) 
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Using the last equations we now prove (2.5) in three steps. First, we give a representation of 
yn,n' ,m i^ terms of a functional of Second, we derive the asyniptotics of yn+i,n+i' ,m 

{td") for 

large n using this representation. Last, we plug this asymptotics into (4.4 1. 

Step 1 (Exact representation of yn,n' ,m(t)). In this step we show that for m < n <n' 



n~m n — m 



-rn — k' 



k=0 k'=0 



(4.6) 



and for m = n < n' 

where S„ n' is Kronecker's S. 



yn.n',m{t) = ^„,„'J/„(t) - yn{t)yn'{t), (4.7) 



Proof. For (4.6), observe that lo„(i) ^^d Y^^i ,_ {t) are independent, given T,„. Moreover, it is 
clear that E[ro„(i)|T„] = y„_,„(c-"(t - r„))lt>T„. Hence, by 

cov[yo„(t),yo„i„,_(i)] 

= E[COV[ro„(i),5^o„i„,_(0|7^m]] 

+ cov[E[yo„ (i) IT™] , E[yo,„i,., W 

1y«-™(c-™(i - T,„))lt>T„, y„'_™(c-™(i - T„))lt>T„J 



n--m n —m 



-rn — k' 



k=0 k'=0 



For (4.7), note that m = n implies that YQ^{t)YQ^i^, _^{t) ~ Sn^n'YQ^{t) and the result follows. □ 

Step 2 (Asymptotics of yn+i-n+i' ,m{tc^) for large n). The aim of this step is to establish that for 
^2, and m € Z+, 



oo oo 



yn+i,n+i' ,m 



(tc")™6LEE 



a-kCLk'S 



2(m+l) _ 1 
..If.' c ± 



k=0 k'=0 



C2 - 1 



(4.8) 



Proof. For large n, we have tc" > and so, by (4.6), (4.7) and dominated convergence, 



oo oo 

yn+i,n+i' ,m{tc"') ~ afcOfe'tLOVIe 

fc=0 fc'=0 

oo oo 



^(tc"-T„) 



„-,i-i' + fe' 



(4.9) 



k=0 fc'=0 
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Now, for all £ Z, using (4.5), for n large enough, 

1 



nr 



^— n— _7+/ ^— n— XT ^ j^— n— j+/ Q—n—j'-\-l i ^— 2n— j/— j'+2/ 



_|_ _|_ ^-2n-2j+2; _|_ 2g-2n-i-j'+2; _|_ ^-2n-2i'+2i-j 



c"""^' + ' + C"""^"+' + c"2n-2j+2; _|_ ^-2n-j-j'+2i _j_ ^-2«-2i'+2;-j 



(4.10) 



^-0 

m 

1=0 
-2n-j~i 



, ^2(m+l) _ 2 



- 1 



Plugging the last expression into (4.9) gives (4.8) 



Step 3 (Combining (4.8) and (4.4)). We write immediately, using j — i — k and / — i' 

m— 1 /c,A;'— 



□ 



(4.11) 



Noting that X^fcLo Wk\c'' < oo, we see that for ai^i'{t) given by (2.6) 

qX„+,(te"),X„+,,(te")] ™ a,,.(t)(^) 
2 /2\2"+'*+«' 



2^^2n+^+^' 2-™(c2™ - 1) 



ai.i'(i) ■ < 



2 - c2 V c 
2"ny2*+^' 



m— 1 

c < 

c = V2, 



C2 - 1 



-2" c'^* 



l2(c2-l)(c2-2) 
This finally shows (2.5) and finishes the proof of Theorem [T] 



> V2. 



□ 



5 Proof of Theorem [2] 

The parameter r is only a rescaling of time. Hence, we can safely assume r = 1 in our proof. Let 



y — {Y{t))t>o be the Aldous-Shields model with parameter c and Xn{t) as in (2.1). Defining 



X^+^it):= E 2''+i-"X„(<), 

n=h+l 

it is important to note that 

Xh+i{t) - #{u e T,,+i : : u ^ v,Y4t) = 1}, 
almost surely; see also Remark |2.5| This implies that we can couple y and Z in the sense that 



t>o 
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By Thcorcm[T] for n < h and i — 0,1,2, . 



h h oo oo 



1=0 j=0 j=0 fc=0 



oo oo 

Z^{tc') ^ Xn+,{tc^) = 2-*+'^/.+.(te'') feoo • 2'^ E E 2a,e-^~-''° 



j=l fc=0 



in probability, for all t > 0. Using the last two limits in the definition of L in (3.1 ) gives the result. 

□ 
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